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We theoretically investigate the terahertz dielectric response of a semiconductor slab hosting an infrared 
photoinduced grating. The periodic structure is due to the charge carries photo-excited by the interference of 
two tilted infrared plane waves so that the grating depth and period can be tuned by modifying the beam 
intensities and incidence angles, respectively. In the case where the grating period is much smaller than the 
terahertz wavelength, we numerically evaluate the ordinary and extraordinary component of the effective 
permittivity tensor by resorting to electromagnetic full-wave simulation coupled to the dynamics of charge 
carries excited by infrared radiation. We show that the photoinduced metamaterial optical response can be 
tailored by varying the grating and it ranges from birefringent to hyperbolic to anisotropic negative dielectric 
without resorting to microfabrication. © 2012 Optical Society of America 
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Active metamaterials exhibiting an externally driven 
electromagnetic response are highly desirable for many 
interesting applications. In the terahertz (THz) fre- 
quency range, many researchers have theoretically and 
experimentally considered a tunable response by semi- 
conductor inclusions in the metamaterial composite. In 
fact, the semiconductor permittivity can be modified by 
exploiting different excitation mechanisms encompassing 
photocarrier injection [1-5], application of a bias volt- 
age [6] and/or thermal carrier excitation [7-9]. Moreover, 
since a spatially modulated optical beam can induce in 
a semiconductor either dielectric or metallic response to 
THz radiation [10], a complex effective metallo-dielectric 
structure can be achieved without time-consuming mi- 
crofabrication processes. For example, T. Okada et al. 
in Ref. [10] considered a periodical grating created by 
an optical pulse on a Si prism showing that the photo- 
induced structure displays metallic properties and, as a 
consequence, it is able to support surface waves. 

In this Letter, we suggest a novel class of infrared 
(IR) driven tunable dielectric metamaterials for THz 
radiation. Here, a periodic dielectric grating is photo- 
generated by the interference pattern of two tilted IR 
plane waves with frequency within the semiconductor 
absorption band; therefore, as opposed to standard pho- 
tonic structures, the bulk grating is amenable to be re- 
configured since one can mold the photo-induced carrier 
density by changing the IR illumination. More precisely, 
the grating depth and period can be straightforwardly 
modified by changing the intensity and incidence angle 
of the two driving IR plane waves, respectively. In the 
situation where the grating period is much smaller than 
the THz wavelength the overall structure behaves as an 
homogeneous medium whose effective permittivity ten- 
sor can be evaluated through first principle simulations. 



We show that the effective dielectric response is gener- 
ally that of a uniaxial crystal with optical axis along 
the grating and that, in a specific range of IR radiation 
intensity, its ordinary and extraordinary principal per- 
mittivities have different signs. As consequence the con- 
sidered structure can act as a hyperbolic material. This 
kind of metamaterials show very remarkable properties 
related to the hyperbolic nature of their isofrequency 
contour such as, for example, subwavelength imaging 
through hyperlensing [11]. Some of the possible realiza- 
tions of hyperbolic media in the THz frequency range in- 
clude homogenous naturally-occurring materials (such as 
triglycine sulfate [12]), periodic array of aligned carbon 
nanotubes [13] and multilayer graphene structure [14]. 

The proposed tunable photo-induced dielectric meta- 
material is sketched in Fig. 1. A semiconductor infinite 
slab of thickness L is illuminated by two IR tilted plane 
waves (at a frequency cj( ir )), namely E^ IR ^ (x, z, t) = 



x - u)( m H\e y with E ( LR > = y/I/(2e c). 
Here r = xe x + ze z , / is the radiation intensity, k.\ = 
k<< 1R )(-sm0e x + cos6e z ) and k^ IR) = k^ lR \s'm 0e x + 
cos0e z ) where k^ lR ^ = o/ IR ) /c and is the incidence an- 
gle and (eo and c are the vacuum permittivity and speed 
of light, respectively). The IR beams photo-excite elec- 
trons in conduction band and such excitation depends on 
the overall intensity distribution and, as a consequence, 
plane waves interference spatially modulates the charge 
carrier density N(x, z) along the x-axis. The steady state 
profile of the photo-generated carrier density is given 
by [15] 
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where D is the diffusion coefficient, r(N) is the density- 
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Fig. 1. (Color on-line) Sketch of tunable photo- induced 
THz dielectric metamaterial slab. The two infrared (IR) 
plane waves (with incidence angles ±6) induce a dielec- 
tric grating along the x-axis affecting the THz radiation. 
The terahertz (THz) plane wave normally impinges on 
the semiconductor slab of thickness L. 

dependent electron-hole recombination time, e^ lR \N) 
is the density-dependent semiconductor permittivity at 
the IR frequency and E^^^x^z) is the amplitude of 
the monochromatic IR field within the semiconductor 
slab (h is the reduced Planck constant). The recom- 
bination rate can be expressed by the polynomial ap- 
proximation l/r(N) = A + BN + CN 2 where the first 
term A is the nonradiative recombination rate coeffi- 
cient whereas the terms with B and the C describe 
radiative transitions and Auger recombinations, respec- 
tively [15]. The IR dielectric permittivity e^ IR ^(7V) can 
be described by adopting the microscopic model in Ref. 
[16,17]. This model is obtained in the quasi-equilibrium 
regime through the free-carrier approximation and con- 
sidering some relevant effects such as the Urbach tail and 
the bandgap renormalization. The incident THz field is 
a plane wave (at a frequency o/ THz )) normally imping- 
ing onto the semiconductor slab interface (see Fig. 1), 
i.e. E™ = E™ cos[k(™ z h-w(™ z H] where E™ 
lies in the x-y plane and /c( THz ) = cj( THz )/c. We can 
described the semiconductor THz dielectric response by 
the Drude dielectric model 



c(THz) 



e 2 r/m*N 



e cj(™ z )(l-^( THz )r) 



(2) 



where is the background dielectric constant, e is the 
electron charge unit, m* is the effective mass of electrons 
in semiconductor and r is the free electron relaxation 
time. Note that by externally changing the intensity / 
and the incidence angle of IR illumination, the carrier 
density TV is locally modified (see Eq. (1)) and this in 
turn induces a spatial modulation of e^ THz ^ (see Eq. (2)). 

The IR driven functionalities of the semiconductor 
slab were investigated through 2D full-wave simulations 
performed with the comsol RF module [18]. THz linear 
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Fig. 2. (Color on-line) (a) Transmissivity T± (Tj|) and 
Reflectivity R± (R\\) of an ordinary (extraordinary) 
wave as a function of the incident IR intensity /. (b) Ef- 
fective dielectric permittivity (e||,ej_) retrieved by trans- 
missivity and reflectivity of panel (a) as a function of the 
incident IR intensity /. (c) Magnification of the shadow 
region of panel (b) where the THz response is hyperbolic. 



and IR nonlinear Maxwell equations with e^ THz ^ given 
by Eq. (2) and e^ IR ^ calculated according to [16, 17] 
were coupled with Eq. (1) for the carrier density pro- 
file. The transmissivity and reflectivity coefficients of the 
slab were evaluated, from which its effective dielectric 
permittivity can be retrieved [19]. 

In our simulations, we have set L = 2.5 /im and 
= 0.04 rad, u/ IR ) = 2168 THz, o/ THz ) = 18.85 THz 
and we have used typical parameters of GaAs bulk at 
room temperature i.e. D = 4 • 10 -3 m 2 /s, A = 10 8 
1/s [16], B = 7.2 • 10- 16 m 3 /s, C = 10~ 42 m 6 /s [20] 
and e b = 13.32 + z8.95 • 10~ 3 [21], r = 3.23 • lO" 13 s, 
m* = 0.067mo (where mo is the electron mass) [22]. We 
have considered that two vacuum layers are placed at the 
facets of the GaAs slab for providing IR and THz exter- 
nal illumination and for evaluating transmission and re- 
flection coefficients. We have adopted periodic boundary 
conditions along the x-axis for both THz and IR radi- 
ation and carrier density. Furthermore, at the entrance 
and the exit facets (orthogonal to the z-axis) of the in- 
tegration domain, we have used matched boundary con- 
ditions for THz and IR electromagnetic fields, whereas 
at the entrance and exit facets of the GaAs slab, we 
have required d z N = 0. Note that the integration do- 
main for N coincides with the GaAs slab where photo- 
generation takes place. Here, the period of the induced 
grating is p = ttc/(uj^ sin 0) = 10.86 /im ~ A( THz Vl0 
so that we assume the slab dielectric response to coincide 
with that of a homogeneous medium at THz frequencies. 
Within such a homogenized regime, THz fields with lin- 
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Table 1. Effective THz dielectric behavior of the considered GaAs slab at different IR intensity. 





Birefringent medium 


Hyperbolic medium 


Anisotropic negative dielectric 




< I (KW/cm 2 ) < 1.3 


1.3 < I (KW/cm 2 ) < 1.45 


I (KW/cm 2 ) > 1.45 


Re(e\\) 


13.3-0 


0- -0.4 


< -0.4 


Im(e\\) 


0-6.0 


6.0-6.8 


> 6.8 


Re(e±) 


13.3-0.9 


0.9-0 


< 


Im(e±) 


0-2 


2-2.2 


> 2.2 



ear polarizations parallel (||-polarized) or orthogonal (_L- 
polarized) to the x axis generally have unequal effective 
dielectric permittivities, e|| and e_L, respectively. In Fig. 
2(a) we report the slab transmissivity and reflectivity of 
1 1- and _L- polarized waves as functions of the IR intensity 
/ in the range KW/cm 2 < / < 4 KW/cm 2 . In Fig. 2(b) 
we plot the effective dielectric functions retrieved from 
the complex reflection and transmission coefficients [19] 
whose squared moduli are reported in panel (a). From 
Fig. 2(b) it is evident that we can steer the THz di- 
electric response by varying the intensity of the two IR 
plane waves. Specifically, the considered structure ex- 
hibits three different behaviors, at different values of the 
IR intensity, whose properties are summarized in Table 
1. At low IR intensities (J < 1.3 KW/cm 2 ) the medium 
behaves as an uniaxial medium (i?e(e||) > Re(e±) > 0). 
In a second IR intensity range 1.3 KW/cm 2 < I < 
1.45 KW/cm 2 , the medium birefringence increases to 
the point that Re(e±) > and Re{e\\) < 0, so that 
the effective permittivity tensor is indefinite and the ef- 
fective medium behaves as an hyperbolic medium. The 
IR intensity range where medium hyperbolicity occurs 
is highlighted both in Fig. 2(b) and in Fig. 2(c) (which 
contains a magnified version of panel (b)) with a shad- 
owed stripe. Finally, at higher IR intensities, for / > 
1.45 KW/cm 2 , the structure effectively displays the be- 
havior of an anisotropic negative dielectric metamaterial 
CHe(e_L_) < and Re(e\\) < 0). 

In conclusion, we have shown that a semiconductor 
slab has a THz dielectric effective behavior which can 
externally be tuned through suitable IR illumination, it 
spanning from slight birefringence to anisotropic nega- 
tive dielectric behavior and, in particulary, it is worth 
noting that the semiconductor slab shows a hyperbolic 
response in a specific range of IR intensity. In addition 
to the multidisciplinary value of our results ensuing from 
the mingling of semiconductor and metamaterial physics 
concepts, it is remarkable that our approach can easily 
be generalized to more complex intensity patterns of the 
IR radiation yielding different and more exotic ways of 
achieving averaged effective medium responses. 

This research has been funded by the Italian Ministry 
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